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The interplay between Mott and Anderson routes to localization in disordered interacting systems gives rise to 
different transitions and transport regimes. Here, we investigate the phase diagram at finite temperatures using 
dynamical mean field theory combined with typical medium theory, which is an effective theory of the Mott- 
Anderson metal-insulator transition. We mainly focus on the properties of the coexistence region associated 
with the Mott phase transition. For weak disorder, the coexistence region is found to be similar as in the 
clean case. However, as we increase disorder Anderson localization effects are responsible for shrinking the 
coexistence region and at sufficiently strong disorder (approximately equal to twice the bare bandwidth) it 
drastically narrows, the critical temperature Tc abruptly goes to zero, and we observe a phase transition in the 
absence of a coexistence of the metallic and insulating phases. In this regime, the effects of interaction and 
disorder are found to be of comparable importance for charge localization. 

PACS numbers: 71.27.+a, 71.30.+h, 71.55.-i, 71.I0.Hf, 71.10.Fd 


I. INTRODUCTION 

Mott mechanism of localizatiorP is an emergent phe¬ 
nomenon in which a large local Coulomb repulsion suppresses 
double occupation, which prevents charge transport in a half- 
filled system. Strongly correlated electron materials, such as 
transition metal oxide J® and some organic salts,®®! exhibit 
a Mott metal-insulator transition due to the effectively strong 
Coulomb repulsion that exists between electrons occupying a 
narrow valence band. Below the critical temperature Tc, this 
transition is of first-order and one observes a region where 
metal and insulator coexist.®® 

The presence of disorder also leads to localization of 
electron wave functi ons, a phenomenon known as Ander¬ 
son localization.®®! in this case the energetic mismatch be¬ 
tween neighboring sites prevents charge transport in the lat¬ 
tice. These two mechanisms of localization - Mott and An¬ 
derson - combine in nontrivial ways, sometimes reducing, 
sometimes enhancing each other’s effects. Recently, the in¬ 
terplay between interaction and disorder has received much 
attention, mainly through three different perspectives. First, 
due to the investigation of the many-body localization,®! a 
novel paradigm arose for understanding localization in dis¬ 
ordered and interacting quantum systems at non-zero temper¬ 
ature. Second, very recently, models of disordered and inter¬ 
acting syst ems have been simulated with cold atoms in opti¬ 
cal lattices.®®! Finally, the disorder and the effective interac¬ 
tion strength can be syst emat ically tuned by doping,®!!®® 
or even X-ray irradiation.®®! 

Over the last few decades considerable progress in the de¬ 
scription of strongly correlated materials and Mott metal- 
insulator transition (MIT) has been achieved through dynam¬ 
ical mean field theory (DMFT).®! In this method, a lattice 
model of interacting electrons is mapped to the Anderson 
impurity model with a conduction bath which needs to be 
calculated self-consistently. To describe disorder, the sim¬ 
plest treatment is within the coherent potential approxima¬ 


tion (CPA).®! The CPA can be easily combined with the 
DMFT,®!®! by considering an ensemble of impurities sur¬ 
rounded by an average bath, which is the same for each elec¬ 
tron. This approach thus does not describe the spatial fluctua¬ 
tions associated with the Anderson localization. Near the An¬ 
derson transition the distribution of the local d ensity of states 
(DOS) changes from Gaussian to log-normal, implying 
that its arithmetic average value does not provide a proper de¬ 
scription of the system. The typical medium theory (TMTj® 
provides a simple method which is able to effectively describe 
the Anderson localization. The central quantity in TMT is the 
typical density of states, defined as the geometric average of 
the local DOS,® which plays the role of the order parameter 
for the Anderson localization. The TMT me thod was carefully 
tested for the noninteracting system,®®IHland it was success¬ 
fully applied to the interacting case within the TMT-DMFT 
approach,® elucidating the full nonmagnetic phase diagram 
for the disordered half-filled Hubbard model and the precise 
nature of the Mott-Anderson critical point.®lThe TMT-DMFT 
approach also allows for a spin-dependence analysis of the 
DOS, which enables one to include effects of long-range mag¬ 
netic order in disordered and interacting systems.® 

In this paper, we perform the first TMT-DMFT calcula¬ 
tion at finite temperatures. We explore the entire nonmagnetic 
phase diagram with a particular focus on the effects of disor¬ 
der on the Mott metal-insulator coexistence region. We care¬ 
fully compare the TMT-DMFT and CPA-DMFT results with 
the goal of precisely determining the Anderson localization 
effects, described only within the former method. We find that 
the TMT-DMFT coexistence region is at comparatively lower 
values of the interaction U, while the critical temperature 
is higher than in CPA-DMFT. The width of the coexistence 
region, however, quickly decreases with disorder. At disor¬ 
der strength W ~ 2B, where B is the bandwidth in the clean 
noninteracting system, TMT-DMFT predicts Tc to abruptly 
go to zero, as opposed to the CPA-DMFT solution where the 
coexistence region asymptotically shrinks to a single point as 
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disorder is increased to infinity. In the regime W > 2B the 
MIT takes place at 17 « W, which makes Anderson and Mott 
mechanisms to become equally important for the properties of 
the system. 

The paper is organized as follows. In section [I^ we briefly 
present the TMT-DMFT method for the solution of the dis¬ 
ordered Hubbard model, and the (U, W) phase diagram is 
shown is section uni Sections |IV] and |V] show details of the 
metal-insulator transition in the presence of weak, moderate, 
and strong disorder. Section|Vl]contains conclusions. 


II. TMT-DMFT METHOD 



We consider the Hubbard model with random site energies, 
given by the Hamiltonian 

H =-t (c|^Cj-a.-|-H.c.) + f/y^ni^n4 + y^(£i-/x)nig., 

{ij)<y i 


where cj^ {cia) creates (destroys) an electron with spin a on 
site i, riia = c\^Cia, t is the hopping amplitude for near¬ 
est neighbor sites, U is the on-site repulsion and £i is the 
random on-site energy, which follows a uniform distribution 
P{£) of width W, centered in = 0. We study the half-filled 
particle-hole symmetric lattice by setting the chemical poten¬ 
tial /r equal to 17/2. In general, transition metal oxides and or¬ 
ganic salts described by the Hubbard model can exhibit both 
antiferromagnetic and paramagnetic Mott insulating phases. 
In this work, we focus on the paramagnetic solution which is 
present even at zero temperature in frustrated lattices. 

Within TMT-DMFT, the lattice model describing a dis¬ 
ordered correlated system is mapped onto an ensemble of 
single-impurity problems, corresponding to sites with dif¬ 
ferent values of the on-site energy, each being embedded 
in a typical effective medium that needs to be calculated 
self-consistently. The TMT-D MFT self-consistent procedure 
can be summarized as followsiSHH By considering an initial 
guess for the (typical) bath A(a;) surrounding the impuri¬ 
ties, we solve an ensemble of impurity problems, which give 
us local Green’s functions G{lo, £i) from which local spectra 
p{uj,£i) = —^lmG{uj,£i) are obtained. The typical DOS is 
then calculated by the geometric average of the local spectra. 


Ptyp(w) = exp 


d£P{£)\np{uj,£) , 


and the typical Green’s function is obtained through the 
Hilbert transform, Gtyp(a;) = doj ^ . For lattices 

with semicircular DOS, po{uj) = in the 

clean non-interacting limit (Bethe lattice with infinite coor¬ 
dination number), the self-consistent loop is closed by calcu¬ 
lating a new bath according to A(a;) = t^Gtyp{oj). To solve 
the single-impurity problems, in this wor k we u se the iterative 
perturbation theory (IPT) on the real axis.l^^® In this case we 
do not need analytic continuation. This is an important advan¬ 
tage of this method since the TMT self-consistency relation is 
based on the local DOS. 


FIG. 1. (Color online) {U, W) phase diagram obtained within TMT- 
DMFT for the disordered Hubbard model at T = 0.008. The de¬ 
scription of the different symbols/colors used is given in the text. 


III. PHASE DIAGRAM 


Fig. presents the TMT-DMFT phase diagram of the dis¬ 
ordered Hubbard model obtained at a small temperature, T = 
0.008. Here and throughout the paper we define the non¬ 
interacting bandwidth B = 4f as the unit of energy. In 
the phase diagram, the black and pink circles correspond 
to the metallic and the insulating spinodal lines of the first- 
order Mott phase transition; these two lines delimit the metal- 
insulator coexistence region. The green triangles indicate a 
transition between a metal and a Mott insulator in the absence 
of a well defined coexistence region (see Section |V]for de¬ 
tails), while the blue stars indicate a transition between a metal 
and a correlated Anderson insulator. Finally, the red squares 
correspond to a crossover between the two insulators, which 
takes place at FF « (7. 

To differentiate the phases and build the phase diagram, we 
have analyzed the behavior of the typical DOS at the Fermi 
level (ptyp(O)), the frequency integrated typical DOS (N) and 
the site occupation as a function of the on-site energy. As an 
example, these quantities are presented in Fig. for the par¬ 
ticular case of 17 = 1.75 and T = 0.008. For this set of 
parameters, as disorder W increases, the system goes from 
the Mott insulator to the Anderson insulator, crossing an in¬ 
termediate metallic phase (see Ref.|40l for example, for a dis¬ 
cussion about the presence of an intermediate metallic phase 
when disorder increases). The Mott insulator is characterized 
by a gap in the typical DOS iptyp{0) = 0) and a finite fre¬ 
quency integrated typical-DOS N (see panel a), as well as a 
single occupation of all sites (see panel b). The metallic phase, 
on the other hand, features a quasi-particle peak in the typical 
DOS, nonzero integrated DOS N, and a variable site occu¬ 
pation Tii. Finally, the correlated Anderson insulator shows 
a vanishing typical DOS, indicating that all the states are lo¬ 
calized and as s uch d o not contribute with spectral weight to 
the typical DOS .1^1221 For this reason, the frequency integrated 
typical DOS goes to zero when the system approaches the An¬ 
derson insulator, and thus can be used as an order parameter 
that signalizes this transition. Furthermore, within the TMT- 
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FIG. 2. (Color online) TMT-DMFT results for 17 = 1.75 and 
T = 0.008. According to the phase diagram of Fig. [T] as disor¬ 
der W increases, the system goes from the Mott insulating phase 
to the metallic phase and finally to the Anderson insulating phase. 
These transitions are identified (and the phase diagram is built) by 
looking at the behavior of the quantities shown in the two panels of 
the present figure: (a) frequency integrated typical DOS, N, and typ¬ 
ical DOS at the Fermi level, ptyp(Q), as a function of W, and (b) site 
occupation per spin as a function of the site energy, normalized by 
the disorder distribution width, W. The inset shows an example of 
the typical DOS in the metallic phase (red long-dashed line), as well 
as in the Mott (black dashed line) and the Anderson (green solid line) 
insulating phases. 


DMFT the Anderson insulating phase corresponds to a two- 
fluid phaseJ^it consists of empty and doubly occupied sites, 
characteristic of non-interacting Anderson insulators, as well 
as singly occupied sites, characterizing Mott localized states 
(see the results for FF = 4 in panel b). 

We find good agreement bet ween our diagram and others 
known in the literature at T = O l^^ l ^^ lThe most relevant effects 
of finite, but small temperature are over the Mott coexistence 
region, which spans over a smaller range of U in comparison 
with the r = 0 case. The real axis IPT impurity solver makes 
it possible to solve TMT-DMFT equations for a broad range 
of parameters and several temperatures. In the following, we 
concentrate on the range of parameters near the phase transi¬ 
tion, and, in particular, near the coexistence region of metallic 


and insulating solutions. 


IV. MOTT TRANSITION FOR WEAK AND MODERATE 
DISORDER W <2B 


In this Section we analyze the coexistence region for weak 
and moderate disorder, which corresponds to FF < W*, 
W* « 1.7. At this regime, the critical U for the Mott transi¬ 
tion is greater than the disorder strength. Although the phase 
transition described within TMT-DMFT is qualitatively simi¬ 
lar as that of CPA-DMFT, some Anderson localization effects 
are already observed. 


A. Coexistence region 

To obtain the coexistence region within CPA-DMFT or 
TMT-DMFT, for a fixed temperature T < Tc, we start from 
a metallic initial bath and increase U to find Uc2, which cor¬ 
responds to the interaction value at which p{0) goes to zero, 
indicating the disappearance of the quasi-particle peak in the 
DOS. Alternatively, when starting from an insulator, by de¬ 
creasing U we find Uc\ where p(0) becomes finite, indicating 
the closure of the gap at the Fermi level. This procedure al¬ 
lows us to obtain hysteresis curves of /9(0) as a function of 
U, which enclose the coexistence region (see Fig.[^for exam¬ 
ples of these hysteresis curves). For a given W, we can repeat 
this procedure for different temperatures and determine the 
two spinodal lines, Uci{T) and Uc2{T), defining the coexis¬ 
tence region. The temperature at which the two spinodal lines 
merge gives the critical temperature, T^, which corresponds to 
a second order critical end point. 

Figure shows the coexistence region obtained as de¬ 
scribed above for the clean case (W = 0) and for a disor¬ 
dered system (W = 0.8), both within TMT-DMFT and CPA- 



u u u 


FIG. 3. (Color online) Hysteresis curves for the DOS at the Fermi 
level obtained by increasing and decreasing 17 at a fixed temperature 
T = 0.01. The curves enclose the coexistence region. The open 
squares were obtained within CPA-DMFT, while the filled circles 
correspond to TMT-DMFT results. 
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FIG. 4. (Color online) Spinodal lines enclosing the coexistence 
region for the clean system (W = 0) and the disordered case 
{W = 0.8) obtained within both TMT-DMFT and CPA-DMFT. 



FIG. 5. (Color online) Average DOS obtained within CPA-DMFT for 
different values of disorder at fixed temperature T = 0.01. Disorder 
broadens the bands in both the metallic (left panel, for (7 = 1) and 
the insulating (right panel, for (7 = 3) phase. 


DMFT. According to our results, when disorder is added to the 
system the critical U at which the transition occurs increases 
in comparison with the clean case. This happens because the 
general effect of disorder is to broaden the bands, as shown 
in Fig. when the CPA-DMFT calculation is performed both 
inside the metallic and the insulating phase. Another general 
effect of disorder seen in the results of Fig. [^is that the tem¬ 
perature of the second order critical point decreases with dis¬ 
order, in agreement with previous CPA-DMFT calculations.!^ 
These general consequences of disorder do not depend on the 
inclusion of Anderson localization effects, since they are ob¬ 
served even within CPA-DMFT approach. 

To carefully study the effects of Anderson localization we 
compare the results obtained within TMT-DMFT with those 
of CPA-DMFT. As can be seen in Fig.|^for W = 0.8, the crit¬ 
ical U at which the transition occurs is smaller within TMT- 
DMFT than within CPA-DMFT. Moreover, a narrower coex¬ 
istence region is observed within the former. To understand 
these results, one should consider that the wave function local¬ 
ization starts at the band edges and that localized states do not 
contribute with spectral weight to the typical DOS. For these 
reasons, in the presence of Anderson localization narrower 
bands are observed in comparison with CPA-DMFT results, 
both in the metallic and the insulating phase, as can be seen 



FIG. 6. (Color online) Typical (TMT-DMFT) and average (CPA- 
DMFT) values of the DOS as a function of frequency showing that 
Anderson localization effects start at the band edges, since both lo¬ 
calized and extended states contribute to the average DOS, while 
only extended ones contribute to the typical DOS. Left panel shows 
results for U = 1.5, while those in the right panel are for U = 1.6, 
both at T = 0.01. 


in Fig. 1^ This is the opposite effect as that described in the 
previous paragraph regarding the effects of adding disorder 
to a clean system. As a consequence, the coexistence region 
within TMT-DMFT is seen in between that of a clean system 
and that obtained within CPA-DMFT for the same value of 
disorder. 

From the TMT-DMFT and the CPA-DMFT hysteresis 
curves shown in Fig. we see that the Anderson localization 
effects over the coexistence region become more important 
as the disorder increases. As W approaches W* « 1.7, the 
width of the TMT-DMFT coexistence region vanishes and we 
were not able to observe the hysteresis even at the lowest tem¬ 
peratures T = 0.005 (see section [V|). In contrast, in the CPA- 
DMFT solution,!^ the coexistence region with finite small Tc 
is observed even for very large W. 


B. Crossover regime and the critical temperature Tc 

As seen in Fig.|^ the coexistence region shrinks as disorder 
increases, making it difficult to obtain the critical temperature 
Tc from the merging of the two spinodal lines. One alternative 
is to determine Tc from the results obtained above it, that is, 
in the crossover region between metal and insulator. This was 
shown to be possible in the clean case and in the present work 
we extend this analysis to the disordered system. 

The quantum Widom line (QWL) associated with the Mott 
transition is defined in Refs. l4TH44l in analogy with the classi¬ 
cal Widom lind^ as the instability (crossover) line above the 
critical end-point {Uc,Tc). It starts at the critical end-point 
and goes to higher temperatures (above the coexistence re¬ 
gion) as a continuation of the first-order phase transition line. 
It is associated with the (zero temperature) quantum critical 
point, which is masked by the coexistence region in the case 
of the Mott transition. The QWL can be defined from the free 
energy functional and can be used to determine 

Tc from the behavior at higher temperatures, as explained (for 
the clean case) in Refs. 1^ and l42l With the objective of ap- 





























5 


plying the QWL analysis to obtain Tc in the disordered case, 
here we review this procedure. 

The Landau free energy functional of the Hubbard model 
as a functional of G{iuJn) is given by 


FL\G(iU}n)\ — —Tt^ G^{iUJn) + Fimp[G{iuJn)], 

n 

where the first term represents the energy needed to form the 
bath around a given site and the second term describes the en¬ 
ergy of the electron at the impurity level surrounded by the 
bath, that is, the free energy of the single-impurity problem. 
The DMFT (TMT-DMFT) equations are obtained by mini¬ 
mizing FL[G{iuJn)] with respect to G{iujn)- 

The curvature A of the above free energy functional with 
respect to U is hnite and minimal along the crossover line 
and is zero at the second order critical point. This curva¬ 
ture can be identified with the convergence rate of the itera¬ 
tive DMFT c alculation,l^il^ that is, A(C/, T) corresponds to the 
slope of the convergence rate ln{Im[G^*‘^(0) — (0)]} 

as a function of the step it of the iterative calculation. Re¬ 
peating the calculation for different values of T, we obtain the 
curve Amin = A(T)|( 7 », where U* is the point at which A 
is minimum for a given T. This line can be extrapolated, to 
Amin|r=Tc = Oj since the curvature of the free energy func¬ 
tional is zero at the second order critical point. 




FIG. 7. (Color online) QWL analysis for the disordered system with 
W = 0.8 described by TMT-DMFT. See the text for the explanation 
of the results in each panel. 


The procedure is illustrated in Fig.j^for the disordered sys¬ 
tem with W = 0.8. For each value of U, we obtain the free 
energy curvature A from the convergence rate of the typical 
Green’s function through the iterative steps, as presented in 
(a) for T — 0.025. For a fixed temperature and different val¬ 
ues of U, we obtain the corresponding A(C/) |t curve. Repeat¬ 
ing this procedure for different temperatures, we obtain the 
set of curves A(?7)|r presented in Fig. |^. The minima Xmin 
of these curves are shown in panel (c), and we obtain Tc as 



U 


FIG. 8. (Color online) QWL and coexistence regions obtained within 
TMT-DMFT (a) and CPA-DMFT (b) for different values of disorder 
(CPA-DMFT coexistence regions for FF > 1.2 were obtained from 
Ref. I26> . The horizontal lines represent Tc obtained from the cor¬ 
responding QWL, calculated as exemplified in Fig.|^(c). The inset 
shows these Tc values as a function of disorder. 


the temperature at which Amm = 0. Finally, panel (d) shows 
the crossover line obtained from data in panel (b), Tc obtained 
through the QWL analysis (gray horizontal line), and the two 
spinodal lines. We conclude that the Tc calculated from the 
QWL analysis coincides with the Tc obtained from the merg¬ 
ing of the two spinodal lines that dehne the coexistence re¬ 
gion. 


In Fig. we show the QWL and the critical temperatures 
obtained from them as we vary the system disorder, both 
within TMT-DMFT and CPA-DMFT. For disorder strengths 
W > 1.6, we find a nonlinear behavior of the TMT-DMFT 
convergence rate as a function of the iteration step; we were 
thus unable to use the QWL analysis discussed to evaluate 
Tc for very large disorder. For W < 1.7, both methods pre¬ 
dict that Tc decreases when W increases (see also the inset in 
Fig. 1^). The critical temperature Tc is higher within TMT- 
DMFT than within CPA-DMFT, although the coexistence re¬ 
gion becomes (very) narrow in the presence of Anderson lo¬ 
calization effects (TMT-DMFT results). However, Tc always 
remains finite within CPA-DMFT even for very large disorder 
strengtfPS whereas we do not observe the coexistence region 
for W > 1.7 in TMT-DMFT (see next section). Our numeri¬ 
cal TMT-DMFT solution indicates that the Tc abruptly drops 
to zero as the coexistence region disappears for W « 1.7. 
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FIG. 9. (Color online) Results obtained within TMT-DMFT for W = 
2.0 at T = 0.01. Panel (a) presents the typical DOS at the Fermi 
level obtained by increasing U (black circles) and decreasing U (red 
stars); no coexistence region is observed. In the same panel we can 
also see the frequency-integrated typical DOS N as a. function of 
U. The inset shows the typical DOS as a function of frequency for 
U = 2.08 (black solid line) and U = 2.10 (red dashed line). Finally, 
panel (b) presents the occupation number per spin as a function of 
the site energy as the transition is approached. 


V. MOTT-ANDERSON TRANSITION FOR STRONG 
DISORDER W>2B 

Within the TMT-DMFT calculation, as we increase disor¬ 
der, the value of the critical U becomes closer to the disorder 
width W. For U W 2B both Mott and Anderson routes 
to localization become equally relevant, and it becomes the 
most difficult to precisely understand the mechanism of the 
MIT. In Fig.we show the results for FF = 2.0 at T = 0.01. 
The transition is seen to take place at (7 « 2.09. Moreover, if 
we look at the results for the typical DOS at the Fermi level 
when U increases, as well as when U decreases (see panel 
(a)), we observe no hysteresis, even if we decrease the tem¬ 
perature down to T = 0.005, in contrast to the results shown 
in Fig. Since ptypiO) becomes zero, the system certainly 
goes through a MIT - but to what type of insulator does the 
system go to? 

To answer the question above, we hrst look at the frequency 


integrated typical DOS N, which can be considered an order 
parameter in the case of the Anderson transition, as discussed 
in the beginning of the paper. As can be seen in Fig. [^a), 
N becomes very small but is still finite when ptyp{0) —1 0, 
suggesting that the transition is not of the Anderson type. The 
nature of the transition can finally be confirmed by analyzing 
the occupation number per spin rii as a function of the site 
energy close to the transition, which can be seen in panel (b). 
As U increases towards the MIT, all sites become singly occu¬ 
pied, which is a characteristic of the Mott insulator. Although 
of the Mott type, the Hubbard subbands are strongly reduced 
for this value of W, as can be seen in the DOS presented in 
the inset, which is consistent with our expectation that both 
Mott and Anderson routes to localization are relevant in this 
regime of t/ « FF. 

Interestingly, our analysis of Fig. [^suggests that for FF = 
2.0 there exists a transition between a metal and a Mott in¬ 
sulator in the absence of a coexistence region. Indeed, ac¬ 
cording to the phase diagram (Fig. [^, the same behavior is 
observed in a small range around [/ « FF « 2. According 
to Figs.j^and TMT-DMFT predicts the coexistence region 
to become (very) narrow when the system is in the U < W 
regime and disorder increases. When the system enters the 
U ^ W regime, the two spinodal lines seem to merge and 
no coexistence is observed, suggesting that Tc abruptly goes 
to zero due to the Anderson localization effects. Our results 
are in general agreement with the T = 0 phase diagram of 
Ref [2] while presenting much more detailed analysis of the 
MIT with the vanishing coexistence region. 

For FF > 2.3, one can find a direct crossover between the 
two insulators, Mott and correlated Anderson, without an in¬ 
termediate metallic phase; this crossover is represented by red 
squares in our diagram of Fig. [T] To distinguish between the 
two insulators, we have looked at the occupation number as a 
function of site energy, as exemplihed in Fig. Our results 
show that when FF < C/ all the sites are single occupied, char¬ 
acterizing a Mott insulator; when FF > (7, on the other hand, 
there are sites with energy larger than (7/2, which are empty, 
sites with energy smaller than —(7/2, which have double oc¬ 
cupancy, and also sites occupied with one electron, character¬ 
izing the two fluid behavior of the correlated Anderson insu¬ 
lator. According to these results, as might have been expected 
from the two fluid picture of the Mott-Anderson insulatopISl 
the crossover between the two insulators is seen to take place 
at WKiU. 


VI. CONCLUSIONS 

In this work, we studied Mott and Anderson routes to local¬ 
ization by using a combination of dynamical mean field the¬ 
ory (DMFT) and typical medium theory (TMT) to solve the 
disordered Hubbard model. According to our TMT-DMFT 
results, Anderson localization has important effects near the 
Mott transition, specially on the coexistence region of metal¬ 
lic and insulating phase that exists below a critical temperature 
Tc. In the presence of small and moderate disorder FF, the 
TMT-DMFT transition is qualitatively similar as in the CPA- 







7 


DMFT case (which does not describe localization due to dis¬ 
order), and the main precursors of the Anderson localization 
are seen in the narrowing of the coexistence region in com¬ 
parison with CPA-DMFT. As the disorder further increases, 
for W > 2B (where B is the bandwidth for (7 = FF = 0), 
the transition occurs at [/ « VF and our results indicate that 
Anderson and Mott routes to localization become equally im¬ 
portant. The critical temperature Tc abruptly goes to zero for 
FF = FF* « 1.7B. For 1.7B <W^U< 2.W the typi¬ 
cal DOS at the metal-insulator transition is strongly reduced, 
but the states are nearly half-filled irrespective of the on-site 
energy, indicating dominantly Mott character of the MIT, al¬ 
though no coexistence region is observed. For even larger dis¬ 
order, FF > 2.3, there is a crossover between the Mott and the 
correlated Anderson insulator. 

The observation of a Mott transition without a coexistence 
region suggests that the nature of the transition has changed 
from first to second order as disorder increases. For the clean 
system, it has been showiJSIthat at T just above Tc the resis¬ 
tivity as a function of temperature shows a scaling behavior. 


which is compatible with an assumption of quantum critical¬ 
ity. In other words, despite the presence of a coexistence re¬ 
gion between the metallic and the Mott insulating phase at 
small temperatures, at intermediate temperatures the system 
seems to be controlled by a hidden quantum critical point. 
Very recently an experimental work on K-organics has con¬ 
firmed the presence of this quantum critical regime at interme¬ 
diate temperatures.® In this respect, it will be very important 
to compare the TMT-DMFT phase diagram and charge trans¬ 
port with the experiments on disordered correlated systems. 
Preliminary results,® on introducing disorder by X-ray irradi¬ 
ation, show that Uc indeed increases with disorder, while Tc 
also decreases, and seems to vanish at some finite disorder. 
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